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PREFACE 

The  US  Air  Force  Environmental  Technical  Applications 
Center  (USAFETAC)  prepared  this  report  in  answer  to  a 
request  from  the  Argonne  National  Laboratory,  Argonne, 
Illinois.  The  information  is  provided  in  support  of  the 
Air  Force  Weapons  Laboratory  (AFWL)  Project  19008W15, 
“Environmental  Implications  of  Airborne  Hydrogen  Fluoride 
(HF)  Laser  Operations . “ 

Argonne  National  Laboratory  requested  USAFETAC  support 
in  the  belief  that  USAFETAC  maintains  a working  computer 
diffusion  model  for  gases  released  from  elevated  line 
sources.  USAFETAC  does  not  have  a working  computer  model; 
therefore,  this  report  presents  generalized  conclusions 
based  upon  a thorough  search  of  references  that  present 
information  relevant  to  this  problem. 

In  the  event  that  this  report  is  incorporated  into 
another  report  by  the  requester  or  any  other  agency, 
request  that  USAFETAC  be  furnished  a copy  of  tae  new 
report  in  all  cases  where  such  dissemination  is  not 
prohibited. 

USAFETAC  prepared  this  report  for  a specific  purpose; 
therefore,  any  further  application  of  this  information 
should  be  undertaken  with  caution.  Work  on  this  report 
progressed  under  rigid  constraints  of  resources  including 
time,  personnel,  equipment,  and  data.  Department  of 
Defense  agencies  and  their  contractors  should  contact 
USAFETAC  directly  for  aid  in  assessing  the  applicability 
of  this  material  for  their  purposes.  Other  prospective 
users  should  contact  professional  environmental  analyrts 
in  the  National  Oceanic  and  Atmospheric  Administration 
(NOAA)  or  private  industry  for  similar  assessment  service. 
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AN  AIR  POLLUTION  ASSESSMENT  OF  HYDROGEN  FLUORIDE 
Introduction 

The  purpose  of  this  report  is  to  estimate  the 
diffusion  characteristics  and  downwind  concentrations 
of  hydrogen  fluoride  (HF)  after  it  is  released  from  an 
airplane.  The  author  conducted  an  exhaustive  search  of 
USAFETAC's  in-house  technical  library  to  determine  a 
method  or  methods  adaptable  to  this  problem. 

Several  authors  have  modified  Sutton's  basic  diffu- 
sion equations  for  microscale  analyses.  Applications 
from  three  sources  are  included  in  this  report.  The 
most  useful  model,  H.  H.  Cramer  (3),  offers  a generalized 
concentration  model  for  a finely  divided  particulate  or 
gas  using  the  Gaussian  distribution  functions.  The 
second  model,  applied  by  D.  B.  Turner  (9) , rewrites 
Sutton's  (1932)  concentration  equation  into  a simplified 
finite  line  source  model  including  an  edge  effects 
term  (EFT) . This  term  is  also  applied  as  the  EET  to 
Cramer's  Generalized  Model  in  this  report.  The  third 
model  is  a reapplication  of  an  equation  by  Smith  and  Hay 
(1961)  using  vertical  turbulence  intensity  (7) . The  US 
Air  Force  empirically  tested  the  diffusion  characteristics 
of  airborne  substances  near  Cedar  Hill,  Texas  and  then 
appropriately  modified  the  original  equation.  Under 
steady-state  conditions,  estimates  can  be  made  using  these 
equations  for  distances  of  up  to  100  kilometers.  These 
line-source  equations  are  briefly  described  in  this  report. 
However,  the  classical  diffusion  equations  are  not  appro- 
priate for  unsteady-state  diffusion  estimates  (5:17,19). 
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Assumptions 

The  state  of  the  art  of  diffusion  estimation  and  the 
limited  amount  of  input  data  require  that  several  assump- 
tions be  made.  They  are: 

(1)  Diffusion  in  the  alongwind  or  x-direction  can 
be  neglected  when  compared  with  a strong  transport  wind 
(5:13) . 

(2)  In  the  vertical  direction  HF  assumes  a statis- 
tical Gaussian  distribution. 

(3)  The  emission  is  an  instantaneous  line  source. 

(4)  Homogeneous  steady-state  conditions  exist  (i.e., 
no  space  or  time  changes  in  wind  or  turbulence) . 

(5)  HF  reflects  perfectly  at  the  surface  (i.e.,  no 
ground  absorption)  and  at  the  height  of  the  mixing  layer. 

(6)  HF  remains  a gas. 

(7)  No  HF  coalesces  with  water  vapor,  or  washes 
or  rains  out. 

(8)  HF  has  approximately  the  same  molecular  weight 
as  air  (i.e.,  no  thermal  buoyancy  or  settling  velocity). 

(9)  Atmospheric  stability  is  neutral  at  all  points 
downwind. 


(10)  Mixing  height  equals  height  of  release. 


tr/ki  ^ 
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(11)  The  mean  wind  direction  is  normal  to  the  air- 
plane's flight  path. 

(12)  No  vertical  wind  shear  at  any  point. 

(13)  The  emission  rate  is  constant. 

Data  Provided 

We  are  given  or  can  calculate  several  pieces  of 
information. 

Mean  Wind  Speed  (m)  = 20  m/s  (45  mph) 

Sampling  will  take  place  near  the  surface  therefore 

z=2  meters 

Four  runs  will  be  made: 


Run  Speed  of  Aircraft  Altitude  of 
Mach  (m/s)  Release  (H) 


Length  of 
Release  (y) 


Summary 


(230)  3 km  (10,000') 

(760)  3 km 

10  km  (33,000') 
10  km 


3,450  m 
10,350  m 
3,450  m 
10,350  m 


Hydrogen  fluoride  (HF)  is  a very  stable  substance 
in  the  atmosphere.  Anhydrous  HF  or  partially  hydrolyzed 
HF  is  completely  soluable  in  the  presence  of  sufficient 
quantities  of  water  vapor.  Thus,  in  order  to  avoid 
coalesence  and  precipitation  scavaging,  no  water  vapor  is 


assumed  to  be  present.  At  various  temperature  and 

pressure  conditions  present  in  the  atmosphere  HF  may 

become  a liquid,  although  it  is  conveniently  assumed  to 

always  be  a gas  here.  Gaseous  HF  is  slightly  lighter 

than  air  but  for  these  computational  purposes  it  will 

have  the  same  mass.  Thus,  only  mechanical  mixing  is 

considered  in  this  report.  Using  the  given  and  calculated 

date  and  the  Jine  source  dispersion  equations  found  in 

Appendix  A,  we  can  estimate  ground-level  and  release-level 

3 

concentration/emission  rate  values  (x/Q  in  sec/m  ) and 
dosages  at  pertinant  distances  downwind. 

Table  1.  Estimated  Ground  Level  Concentrations  (from 
Cramer  (1) ) 


Concentrations  ( s ec/nr ) 


Distances  Downwind 


60  km 

100  km 

*(x/Q)3 

4.0xl0“5(max) 

3.9xl0~5 

*(X/Q)10 

8.0xl0~7 

4 . OxlO”6 

Dosage^ 

4.0xl0~^ (max) 

3.9xl0*"5 

Dosage10 

2. 8xl0**8 

4.0xl0~6 

♦Subscripts  denote  release 

from  3 km  and  10  km. 

Maximum  values  for  releases  at  10  km  altitude  are 
given  for  100  km  downwind.  The  mathematical  peak  occurs 
beyond  that  distance  but  the  model  may  not  be  accurate 
past  100  km.  If  the  mixing  is  below  the  release  height 
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of  the  HF  gas , virtually  no  ground  contamination  will 
occur. 

Consistant  and  accurate  meteorological  inputs  are 
paramount  £or  the  successful  application  of  any  pollution 
concentration  estimation  method.  Among  the  meteorological 
parameters , the  mixing  depths  downwind  from  the  release 
point  are  the  most  important  factors  in  determining  cloud 
expansion  and  ground  contamination. 

The  assumptions  made  at  the  outset  of  this  report 
make  the  expected  concentrations  extreme  worst  case  figures. 
This  conservatism  means  that  under  virtually  all  meteoro- 
logical conditions,  the  expected  concentrations  will  not 
exceed  those  given. 

Discussion 

Using  Cramer's  equation  (3:21)  we  can  calculate 
center  line  concentrations  and  graph  the  results  for 
X/Q  values  at  the  release  height  (Figure  1)  or  at  the 
ground  level  when  released  from  3 km  (Figure  2) . When  HF 
is  released  from  3 km,  the  maximum  ground  level  concent- 
ration occurs  60  km  downwind  and  the  corresponding  x/Q  is 
is  3.22xl0*"5  sec/m3. 

When  HF  gas  is  released  from  an  altitude  of  10  km, 

the  maximum  concentration  theoretically  occurs  200  km 

-5  3 

downwind  and  is  about  10  sec/m  . However,  Cramer's 

steady-state  equation  is  not  valid  beyond  100  km.  At 

—6  3 

this  distance  x/Q  is  about  5.2x10  sec/m  . 

The  EET  (Figure  A-4)  is  a function  of  the  length 
of  the  spray  line,  y,  which,  when  considering  a constant 
emission  time,  varies  directly  with  aircraft  speed.  Thus, 
at  Mach  0.5,  the  concentration  values  at  either  “Jdge  of 
the  3450  meter  spray  line  falls  off  by  only  4%  while  at 
Mach  1.5,  virtually  no  concentration  loss  can  be  noticed 
at  the  spray  line's  edge. 
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Figures  3 and  4 give  the  vertical  term  (VT)  versus 
the  distance  downwind  for  a release  altitude  of  10  km. 
Similar  calculations  for  using  the  edge  effects  term  can 
be  made  for  the  release  height. 

Figure  5 illustrates  the  estimated  center  line  dosage 
term  versus  distance  downwind.  This  is  interpreted  as  the 
amount  of  HF  that  passes  a point  during  an  entire  spray 
episode.  The  EET  (Figure  A-4)  is  applicable  to  the  dosage 
terms  as  well. 

Figures  6 and  7 show  the  expected  ground  level 
concentration  of  HF  when  released  from  an  altitude  of  3 km 
and  10  km  respectively. 

Other  investigations  to  determine  the  actual  downwind 
concentration  from  elevated  line  sources  have  been  con- 
ducted using  empirical  experiments.  Among  them  was  an 
Air  Force  test  at  the  Cedar  Hill/  Texas  television  tower 
(7:171).  The  Air  Force  released  traces  of  zinc  cadmium 
sulfide  from  a low  flying  aircraft  while  samp?ers  were 
placed  at  regular  intervals  downwind  to  48  km.  The 
experimenters  related  cloud  expansion  to  meteorological 
parameters  including  vertical  and  horizontal  turbulence 
and  wind  velocity.  In  a well  developed  turbulent  layer/ 
estimated  ground  level  concentrations  agreed  well  with 
the  mathematical  model  listed  below. 


Sxu/ZlT 


"For  releases  above  this  turbulent  layer/  ground  exposures 
were  much  more  erratic  than  predicted  by  the  model  (1:171)." 


exp 


-[A] 

i 18X4*  i 


(1) 


When  vertical  turbulence  is  unity  the  released  alti- 
tude concentration  can  be  plotted  for  distances  downwind 
(Figure  8).  Cramer's  VT,  calculated  above,  can  be  used 
together  with  the  above  C/Q  to  estimate  ground  level 
concentrations . 
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Turner  (9:41)  developed  an  equation  to  estimate 
dovmwind  concentrations  at  the  ground  from  a finite  line 
source.  The  wind  must  be  normal  to  the  spray  line.  The 


standard  deviation  in  the  vertical,  o . is  taken  from 

7, 


the  calculations  for  Cramer's  model. 

Figure  9 is  an  illustration  of  the  estimated  dovm- 
wind ground  concentrations  from  Turner's  equation: 


X/Q 


2ir  o u 
z 


exp  [ -1/2 (-1/2 (§  )2J 


(2) 


For  distances  greater  than  100  km  downwind,  in  the  meso- 
scale  under  turbulent  conditions,  or  in  the  macroscale, 
the  classical  diffusion  equations  are  not  appropriate. 
Predictions  of  concentration  distributions  are  more 
accurately  made  with  synoptic  forecasts  of  the  movement 
of  large  air  masses.  USAFETAC  does  not  have  the  capability 
to  make  large  scale  turbulence  predictions  at  this  time. 
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Figure  6.  Estimated  Ground  Leve 
Dosage  (x/Q)  versus  Downwind 
Distance  when  HF  Gas  is  Released 
from  an  Altitude  of  3 km. 
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Figure  9.  Concentration/ 
Emission  Rate,  x/Q#  at  Ground 
versus  Distance  Downwind  from 
Turner  at  Two  Release  Altitudes 
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Appendix  A 

CALCULATIONS  FOR  CRAMER'S  GENERALIZED  CONCENTRATION  MODEL 

Craruer  (3:21)  uses  a generalized  mathematical  predic- 
tion model  containing  five  terms  to  calculate  downwind 
line  source  concentrations.  The  model  is  simply  written: 


Concentration  (x)  = CCT  x VT  x EET  x AT  x DT 


(A-l) 


where 


CCT  = Centerline  Concentration  Term 


iiraZ°x 


(A-2) 


The  emission  rater  Q,  is  unknown,  therefore  when  we  solve 
for  X/Q, 


CCT  = 


7roz°x 


(A- 3) 


where 


= standard  deviation  of  vertical  concentration 
distribution 


= *E<*rz> 


^z“xrz(1-e> 


(A-4) 


where 


o'  = UA/3 


(A-5) 


m«ier 


where 


o»  - 


-rrs:^ 

from  A'“  expected,  value, 

. xo  (using  the  medtan  exp 


■ ftV/ytw . 


vertical  diffusion  coefficient  (3:63) 
references  sources  that  show  this  value 
approaches  unity  for  elevated  releases 
under  neutral  conditions. 


x_  = 


vertical  virtual  distance  - when  the 
standard  deviation  of  the  vertical 
concentration  distribution  is  small, 
this  value 


zR  - x 


Rz 


1/3.3  ■ 0.30 


(A-6) 


xRz  * distance  downwind  at  which  the 

standard  deviation  of  the  vertical 
concentration  distribution  is 
measured 

* 0 (since  we  are  interested  only 
in  the  standard  deviation  at  the 
source) 


thus,  as  an  example,  assume  that  we  are  interested  in  a 

4 

distance  10  km  (10  m)  downwind. 


then 


a « .052(20) 
z 


= 500 


104+0. 3-20(1-1)  1 
TJJU) 


(A-7) 


Figure  A-2  plots  all  oz  for  distances  downwind  to  100  km. 
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i 


f 

*jfe 

>r 

A 


thus, 


a * standard  deviation  of  the  downwind  concentration 
distribution. 

a 1 (using  assumptions  (1)  and  (12)) 


CCT/Q 


3.14xl0**4  sec/m3 


(A-8) 


for  a release  altitude  of  3 km  at  a distance  of  10  km 
downwind. 


VT  = vertical  term 


= exp 


[-1/2  <5j£) 2 J + exp  [ - 1/2  (Sii) 2 J 


“ r 2iH  -H-z  2-1  r 2iH  -H+z  2-i 

I (exp  l‘V2< ^ > J + exP  [ “1/2  ( ^ > j 


[2iH  +H-Z  2 
-1/2  ( § ) 


2*1  r 2iH+H+z  2-j 

) J + exp  -1/2  (—2 ) J (a-9) 


The  vertical  term  refers  to  the  expansion  of  the  gas  in 
the  z-direction  where  for  Run  (1) 

H = height  of  release  * 3000  m 
z = height  of  interest  = 2 m (ground  level) 

Hm  * 3000  ^see  assumption  (1)) 
o = 500  m (at  10  km) 

Z 

for  one  iteration 

VT  - exp  [-l/2(3ijjj{j~2)2  J + exp  [-l/2(^.ggl2-)  j ♦ 


T , ,6000-3000-2, 21, T , ,„  ,6000-3000+2, ' 

exp  1^  -1/2  ( jyg ) J+  exp  -1/2  ( ^ ) 


, F , ,,,,6000+3000-2, 21. T , ,„  ,6000+3000+2, 

+ exp  -1/2  ( 533 > J + exp  -1/2  ( 533 ) 


exp (-18)  + exp (-18)  + exp (-18)  + exp (-18) 
+ exp (-162)  + exp (-162) 


(A-10) 


VT  = 4 exp (-18)  + 2 exp (-162) 


(A— 11) 


» 1.52x10  ° 

Only  one  iteration  is  used  in  the  vertical  term  calculation 
because  succeeding  iterations  become  negligibly  small 
relative  to  the  total  vertical  term.  Figure  A-3  relates  the 
vertical  diffusion  term  with  downwind  distances. 


EET  « J**2  — exp  (-l/2p2)  dp 

px  f2Tt 


, , 2.  « r -(2i-l)^>2 

-exp(-l/2p  ) r -p 

Z*P  .£2  x 

Z (2n-3) 

- n=2  J 


(Turner  (9:41))  (A-12) 


where 


px  * ^1  and  p.,  » y2 
c ff 

y y 


when  the  spray  line  stretches  from  y^  to  y2* 


where 


a'( t)  = (a'(r  n(I-)1/5  (A-14) 

Cl  Cl  Q l 

0 

where 

t = emission  time  = 15  seconds 
To  = reference  time  = 600  seconds 

t 

oft(To)  = standard  deviation  of  the 

angle  measured  over  reference 
time  (to) 

- 15®  = .263  radians  (from 
Cramer  (3:53)) 


K 


' i 


a^(T)  =(.263)  (^|)1/5  = .263(.025)  * 

v-4 
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=1.58x10 


ry 


distance  at  which  crosswind  cloud 
stabilization  occurs  downwind  from 
source - 

40  meters  (from  Dettling  (4:12)) 


10*  meters  (from  example  in  Equation 
A-7) 


crosswind  virtual  distance 


23 


B 


Ri- 


ff 
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l» 
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Is 
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tt 


K 


fc 


ft 


m 


1 


II 


^V  " ’"'  ,'^^'<S^*^; 


-SL-  “ xRy 


aA(x) 


(A-15) 


where 


OyR  = standard  deviation  of  the  cross- 


wind  distribution 


x„„  = 


Ry 


distance  at  which  the  standard 
deviation  of  crosswind  concentra- 
tion is  measured 
0 


a = crosswind  diffusion  coefficient 
= 1 


Ad  » azimuth  wind  direction  shear  between 
ground  level  and  release  level 


as 


Iz  tz2-zl> 


(A-16) 


where 


A6 

Iz 


rate  change  of  wind  direction 
from  surface  to  release  height 
(radians/meters) 


Since  the  HF  is  released  well  above  the  gradient  level , the 
angle  between  the  surface  and  gradient  level  is  assumed  tc 
be  45°  as  estimated  by  Sutton  (6:71). 


.14  radians 
3660  meters 


(A-17) 
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0.47x10  radians/meters 


z2  ■ effective  upper  bound  of  cloud 

* 3000  meters 

z^  * effective  lover  bound  of  cloud 

* 2 meters 

A0 1 » 0.47xl0“4 (3000-2)  (A-18) 

= .14 


thus 

ay  = a. 58x10-^(40)  1,2^04^01, 2 V* 

(A-19) 

= 811  meters  (for  a distance  10  km  downwind) 


Finally,  to  calculate  the  limits  of  the  EET  summation  term 

y1  j * distance  from  centerline  of  concen- 
tration 
* 1725  meters 


then 

^1  # 2 1725  . < . /•  .nt 

pl,2  "5IT  2,13  (A~20) 

+2.13 

EET  » / ^ exp  (-1/2 (2.13) 2) dp  (A-21) 

-2.13 
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m 


* 


-i'Zi  £ 


I 


I 


I 

S' 

» 


* .9668  (from  Burrington  (2:273))  for  a sampling 
height  of  2 m when  HF  is  released  from  a height 
of  3 km. 

Figure  A-4  expresses  EET  relative  to  airplane  speed  at 
release  time. 

AT  - cloud  growth  in  the  along  wind  direction 
® 1 (using  assumptions  (1)  and  (13) 


I 

* 


DT  - depletion  term  or  the  loss  of  material  by  decay 
processes 

■ 1 (using  assumption  (5) ) 


I 


Figure  A-3.  Vertical  Term  of 
Cramer's  Concentration  Equation 
for  a Release  Altitude  of  3 Km. 
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